DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems.

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

Subtle transcriptomic response of Eurasian perch (Perca fluviatilis) associated with Triaenophorus nodulosus plerocercoid infection.

Determining the physiological effects of parasites and characterizing genes involved in host responses to infections are essential to improving our understanding of host-parasite interactions and their ecological and evolutionary consequences. This task, however, is complicated by high diversity and complex life histories of many parasite species. The use of transcriptomics in the context of wild-caught specimens can help ameliorate this by providing both qualitative and quantitative information on gene expression patterns in response to parasites in specific host organs and tissues. Here, we evaluated the physiological impact of the widespread parasite, the pike tapeworm (Triaenophorus nodulosus), on its second intermediate host, the Eurasian perch (Perca fluviatilis). We used an RNAseq approach to analyse gene expression in the liver, the target organ of T. nodulosus plerocercoids, and spleen which is one of the main immune organs in teleost fishes. We compared perch collected from multiple lakes consisting of individuals with (n = 8) and without (n = 6) T. nodulosus plerocercoids in the liver. Results revealed a small number of differentially expressed genes (DEGs, adjusted p-value ≤0.05) in both spleen (n = 22) and liver (n = 10). DEGs in spleen consisted of mostly upregulated immune related genes (e.g., JUN, SIK1, THSB1), while those in the liver were often linked to metabolic functions (e.g., FABP1, CADM4, CDAB). However, Gene Ontology (GO) analysis showed lack of functional enrichment among DEGs. This study demonstrates that Eurasian perch displays a subtle response at a gene expression level to T. nodulosus plerocercoid infection. Given that plerocercoids are low-metabolic activity transmission stages, our results suggest that moderate T. nodulosus plerocercoid infection most likely does not provoke an extensive host immune response and have relatively low physiological costs for the host. Our findings illustrate that not all conspicuous infections have severe effects on host gene regulation.

Spatial and intra-host distribution of myxozoan parasite Tetracapsuloides bryosalmonae among Baltic sea trout (Salmo trutta).

Proliferative kidney disease caused by the myxozoan parasite Tetracapsuloides bryosalmonae has been actively studied in juvenile salmonids for decades. However, very little is known about parasite prevalence and its geographical and intra-host distribution at older life stages. We screened T. bryosalmonae among adult sea trout (Salmo trutta) (n = 295) collected along the Estonian Baltic Sea coastline together with juvenile trout from 33 coastal rivers (n = 1752) to assess spatial infection patterns of the adult and juvenile fish. The parasite was detected among 38.6% of adult sea trout with the prevalence increasing from west to east, and south to north, along the coastline. A similar pattern was observed in juvenile trout. Infected sea trout were also older than uninfected fish and the parasite was detected in sea trout up to the age of 6 years. Analysis of intra-host distribution of the parasite and strontium to calcium ratios from the otoliths revealed that (re)infection through freshwater migration may occur among adult sea trout. The results of this study indicate that T. bryosalmonae can persist in a brackish water environment for several years and that returning sea trout spawners most likely contribute to the parasite life cycle by transmitting infective spores.

Integrating population genetics and species distribution modelling to guide conservation of the noble crayfish, Astacus astacus, in Croatia.

The noble crayfish, Astacus astacus, is an indigenous European freshwater species. Its populations show significant declines caused by anthropogenic pressure on its habitats, climate change and the spread of invasive species. Diminishing populations' trends and loss of genetic diversity highlight the need for effective conservation that will ensure their long-term survival. We combined population genetics and species distribution modelling (SDM) to reveal the impact of climate change and invasive species on the noble crayfish, and to guide future conservation programs of current populations. Our study showed that Croatian populations of A. astacus harbour an important part of species genetic diversity and represent significant genetic reservoir at the European level. The SDM results predicted substantial reductions of suitable habitats for A. astacus by the 2070; only 13% of its current potential distribution is projected to remain stable under pessimistic Representative Concentration Pathway (RCP 8.5) emission scenario. Moreover, most of the populations with high genetic diversity are located in the areas predicted to become unsuitable, and consequently have a high probability of being lost in the future. Further, SDM results also indicated considerable decrease of future habitat suitability for invasive crayfish species in Croatia, suggesting that climate change poses a major threat to already endangered A. astacus. The obtained results help in the identification of populations and areas with the highest conservation value which should be given the highest priority for protection. In order to preserve present diversity in areas that are predicted as suitable, we propose assisted migration and repopulation approaches, for enhancing populations' size and saving maximum genetic variability. The result of our research emphasizes once again the benefits of multidisciplinary approach in the modern biodiversity conservation.

Whole-genome sequencing illuminates multifaceted targets of selection to humic substances in Eurasian perch.

Extreme environments are inhospitable to the majority of species, but some organisms are able to survive in such hostile conditions due to evolutionary adaptations. For example, modern bony fishes have colonized various aquatic environments, including perpetually dark, hypoxic, hypersaline and toxic habitats. Eurasian perch (Perca fluviatilis) is among the few fish species of northern latitudes that is able to live in very acidic humic lakes. Such lakes represent almost "nocturnal" environments; they contain high levels of dissolved organic matter, which in addition to creating a challenging visual environment, also affects a large number of other habitat parameters and biotic interactions. To reveal the genomic targets of humic-associated selection, we performed whole-genome sequencing of perch originating from 16 humic and 16 clear-water lakes in northern Europe. We identified over 800,000 single nucleotide polymorphisms, of which >10,000 were identified as potential candidates under selection (associated with >3000 genes) using multiple outlier approaches. Our findings suggest that adaptation to the humic environment may involve hundreds of regions scattered across the genome. Putative signals of adaptation were detected in genes and gene families with diverse functions, including organism development and ion transportation. The observed excess of variants under selection in regulatory regions highlights the importance of adaptive evolution via regulatory elements, rather than via protein sequence modification. Our study demonstrates the power of whole-genome analysis to illuminate the multifaceted nature of humic adaptation and provides the foundation for further investigation of causal mutations underlying phenotypic traits of ecological and evolutionary importance.

Genetic diversity and structure of the noble crayfish populations in the Balkan Peninsula revealed by mitochondrial and microsatellite DNA markers.

BACKGROUND: The noble crayfish (Astacus astacus) is a native European species in decline, with a contracting range and diminishing populations and abundance. Previous studies revealed this species significant genetic diversity in the south-eastern Europe, with populations from the western and the southern part of the Balkan Peninsula being the most divergent. However, sampling of populations from the western part of the Balkans was limited and insufficient for investigating genetic diversity and population divergence for the purpose of conservation planning and management. Thus, the major aim of this study was to fill in this knowledge gap by studying mitochondrial and microsatellite DNA diversity, using 413 noble crayfish from 18 populations from waterbodies in the western part of the Balkan Peninsula. METHODS: Phylogenetic analysis of studied populations and their mitochondrial diversity were studied using COI and 16S sequences and population genetic structure was described using 15 microsatellite loci. RESULTS: Phylogeographic analysis revealed new divergent mitochondrial haplotypes for the populations in the westernmost part of the Balkan Peninsula in the tributaries of the Sava and Drava rivers. Microsatellite data indicated that these populations harbour an important component of genetic diversity within A. astacus. The results suggest that the western part of the Balkans played an important role as microrefugia during the Pleistocene climate fluctuations, allowing the long term persistence of A. astacus populations in this region. These results will also be important to supporting conservation decision making and planning.

Climate change-driven disease in sympatric hosts: Temporal dynamics of parasite burden and proliferative kidney disease in wild brown trout and Atlantic salmon.

Global climate change is altering the abundance and spread of various parasites, which has important consequences not only for host-parasite interactions but also for the relationships between different host species. Here, we focus on the myxozoan endoparasite Tetracapsuloides bryosalmonae that causes temperature-dependent proliferative kidney disease (PKD) in salmonids. We characterized the temporal changes in the parasite load and the severity of PKD signs (renal hyperplasia, haematocrit) in two sympatric populations of wild brown trout (Salmo trutta) and Atlantic salmon (Salmo salar). We found that both the parasite load and disease signs vary considerably between individuals, species, rivers and sampling periods. We showed that Atlantic salmon was able to slow down the initial parasite proliferation rate and subsequently tolerate high parasite burden without obvious disease signs. In contrast, the initial parasite proliferation rate was much higher in brown trout, which was followed by the development of severe PKD signs. Thus, the speed of parasite proliferation, rather than the absolute number of the parasites in the host kidney, may play an important role in interspecific variation in PKD susceptibility. To conclude, this study illustrates the usefulness of temporal perspective for understanding host defence mechanisms and climate change-mediated impacts in the wild.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids.

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

The strength and form of natural selection on transcript abundance in the wild.

Gene transcription variation is known to contribute to disease susceptibility and adaptation, but we currently know very little about how contemporary natural selection shapes transcript abundance. Here, we propose a novel analytical framework to quantify the strength and form of ongoing natural selection at the transcriptome level in a wild vertebrate. We estimated selection on transcript abundance in a cohort of a wild salmonid fish (Salmo trutta) affected by an extracellular myxozoan parasite (Tetracapsuloides bryosalmonae) through mark-recapture field sampling and the integration of RNA-sequencing with classical regression-based selection analysis. We show, based on fin transcriptomes of the host, that infection by the parasite and subsequent host survival is linked to upregulation of mitotic cell cycle process. We also detect a widespread signal of disruptive selection on transcripts linked to host immune defence, host-pathogen interactions, cellular repair and maintenance. Our results provide insights into how selection can be measured at the transcriptome level to dissect the molecular mechanisms of contemporary evolution driven by climate change and emerging anthropogenic threats. We anticipate that the approach described here will enable critical information on the molecular processes and targets of natural selection to be obtained in real time.

Draft Genome Assembly of the Freshwater Apex Predator Wels Catfish (Silurus glanis) Using Linked-Read Sequencing.

The wels catfish (Silurus glanis) is one of the largest freshwater fish species in the world. This top predator plays a key role in ecosystem stability, and represents an iconic trophy-fish for recreational fishermen. S. glanis is also a highly valued species for its high-quality boneless flesh, and has been cultivated for over 100 years in Eastern and Central Europe. The interest in rearing S. glanis continues to grow; the aquaculture production of this species has almost doubled during the last decade. However, despite its high ecological, cultural and economic importance, the available genomic resources for S. glanis are very limited. To fulfill this gap we report a de novo assembly and annotation of the whole genome sequence of a female S. glanis The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a highly continuous draft genome of S. glanis: ∼0.8Gb assembly (scaffold N50 = 3.2 Mb; longest individual scaffold = 13.9 Mb; BUSCO completeness = 84.2%), which included 313.3 Mb of putative repeated sequences. In total, 21,316 protein-coding genes were predicted, of which 96% were annotated functionally from either sequence homology or protein signature searches. The highly continuous genome assembly will be an invaluable resource for aquaculture genomics, genetics, conservation, and breeding research of S. glanis.

Humic-acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding.

BACKGROUND: Next generation sequencing (NGS) technologies are extensively used to dissect the molecular mechanisms of host-parasite interactions in human pathogens. However, ecological studies have yet to fully exploit the power of NGS as a rich source for formulating and testing new hypotheses. METHODS: We studied Eurasian perch (Perca fluviatilis) and its eye parasite (Trematoda, Diplostomidae) communities in 14 lakes that differed in humic content in order to explore host-parasite-environment interactions. We hypothesised that high humic content along with low pH would decrease the abundance of the intermediate hosts (gastropods), thus limiting the occurrence of diplostomid parasites in humic lakes. This hypothesis was initially invoked by whole eye RNA-seq data analysis and subsequently tested using PCR-based detection and a novel targeted metabarcoding approach. RESULTS: Whole eye transcriptome results revealed overexpression of immune-related genes and the presence of eye parasite sequences in RNA-seq data obtained from perch living in clear-water lakes. Both PCR-based and targeted-metabarcoding approach showed that perch from humic lakes were completely free from diplostomid parasites, while the prevalence of eye flukes in clear-water lakes that contain low amounts of humic substances was close to 100%, with the majority of NGS reads assigned to Tylodelphys clavata. CONCLUSIONS: High intraspecific diversity of T. clavata indicates that massively parallel sequencing of naturally pooled samples represents an efficient and powerful strategy for shedding light on cryptic diversity of eye parasites. Our results demonstrate that perch populations in clear-water lakes experience contrasting eye parasite pressure compared to those from humic lakes, which is reflected by prevalent differences in the expression of immune-related genes in the eye. This study highlights the utility of NGS to discover novel host-parasite-environment interactions and provide unprecedented power to characterize the molecular diversity of cryptic parasites.

Highly Continuous Genome Assembly of Eurasian Perch (Perca fluviatilis) Using Linked-Read Sequencing.

The Eurasian perch (Perca fluviatilis) is the most common fish of the Percidae family and is widely distributed across Eurasia. Perch is a popular target for professional and recreational fisheries, and a promising freshwater aquaculture species in Europe. However, despite its high ecological, economical and societal importance, the available genomic resources for P. fluviatilis are rather limited. In this work, we report de novo assembly and annotation of the whole genome sequence of perch. The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a draft perch genome ∼1.0 Gbp assembly (scaffold N50 = 6.3 Mb; the longest individual scaffold of 29.3 Mb; BUSCO completeness of 88.0%), which included 281.6 Mb of putative repeated sequences. The perch genome assembly presented here, generated from small amount of starting material (0.75 ng) and a single linked-read library, is highly continuous and considerably more complete than the currently available draft of P. fluviatilis genome. A total of 23,397 protein-coding genes were predicted, 23,171 (99%) of which were annotated functionally from either sequence homology or protein signature searches. Linked-read technology enables fast, accurate and cost-effective de novo assembly of large non-model eukaryote genomes. The highly continuous assembly of the Eurasian perch genome presented in this study will be an invaluable resource for a range of genetic, ecological, physiological, ecotoxicological, functional and comparative genomic studies in perch and other fish species of the Percidae family.

Quantitative Genetic Variation in, and Environmental Effects on, Pathogen Resistance and Temperature-Dependent Disease Severity in a Wild Trout.

Health after pathogen contact varies among individuals because of differences in pathogen load (which is limited by resistance) and disease severity in response to pathogen load (which is limited by tolerance). To understand pathogen-induced host evolution, it is critical to know not only the relative contributions of nongenetic and genetic variation to resistance and tolerance but also how they change environmentally. We quantified nongenetic and genetic variation in parasite load and the associated temperature-dependent disease among trout siblings from two rivers. We detected a genetic variance for parasite load 6.6 times as large in the colder river. By contrast, genetic variance for disease traits tended to be larger in the warmer river, where the disease was manifested more severely. The relationships between disease severity and pathogen load (tolerance) exhibited plateaus at low pathogen load and stronger steepening slopes at high pathogen load in the warmer river. Our study demonstrates the environmental influence on disease severity, nongenetic and genetic variance for health-damage-limiting host abilities, and the shape of tolerance curves. Environmental variability is predicted to govern the presence and intensity of selection, change the relative contributions of nongenetic and genetic variance, and therefore hamper evolution toward more resistant and tolerant hosts.

Outlier Loci Detect Intraspecific Biodiversity amongst Spring and Autumn Spawning Herring across Local Scales.

Herring, Clupea harengus, is one of the ecologically and commercially most important species in European northern seas, where two distinct ecotypes have been described based on spawning time; spring and autumn. To date, it is unknown if these spring and autumn spawning herring constitute genetically distinct units. We assessed levels of genetic divergence between spring and autumn spawning herring in the Baltic Sea using two types of DNA markers, microsatellites and Single Nucleotide Polymorphisms, and compared the results with data for autumn spawning North Sea herring. Temporally replicated analyses reveal clear genetic differences between ecotypes and hence support reproductive isolation. Loci showing non-neutral behaviour, so-called outlier loci, show convergence between autumn spawning herring from demographically disjoint populations, potentially reflecting selective processes associated with autumn spawning ecotypes. The abundance and exploitation of the two ecotypes have varied strongly over space and time in the Baltic Sea, where autumn spawners have faced strong depression for decades. The results therefore have practical implications by highlighting the need for specific management of these co-occurring ecotypes to meet requirements for sustainable exploitation and ensure optimal livelihood for coastal communities.

Less is more: extreme genome complexity reduction with ddRAD using Ion Torrent semiconductor technology.

Massively parallel sequencing a small proportion of the whole genome at high coverage enables answering a wide range of questions from molecular evolution and evolutionary biology to animal and plant breeding and forensics. In this study, we describe the development of restriction-site associated DNA (RAD) sequencing approach for Ion Torrent PGM platform. Our protocol results in extreme genome complexity reduction using two rare-cutting restriction enzymes and strict size selection of the library allowing sequencing of a relatively small number of genomic fragments with high sequencing depth. We applied this approach to a common freshwater fish species, the Eurasian perch (Perca fluviatilis L.), and generated over 2.2 MB of novel sequence data consisting of ~17,000 contigs, identified 1259 single nucleotide polymorphisms (SNPs). We also estimated genetic differentiation between the DNA pools from freshwater (Lake Peipus) and brackish water (the Baltic Sea) populations and identified SNPs with the strongest signal of differentiation that could be used for robust individual assignment in the future. This work represents an important step towards developing genomic resources and genetic tools for the Eurasian perch. We expect that our ddRAD sequencing protocol for semiconductor sequencing technology will be useful alternative for currently available RAD protocols.

Wild Estonian and Russian sea trout (Salmo trutta) in Finnish coastal sea trout catches: results of genetic mixed-stock analysis.

For responsible fisheries management of threatened species, it is essential to know the composition of catches and the extent to which fisheries exploit weak wild populations. The threatened Estonian, Finnish and Russian sea trout populations in the Gulf of Finland are targets of mixed-stock fisheries. The fish may originate from rivers with varying production capacities, from different countries, and they may also have either a wild or hatchery origin. In order to resolve the composition of Finnish coastal sea trout catches, we created a standardized baseline dataset of 15 DNA microsatellite loci for 59 sea trout populations around the Gulf of Finland and tested its resolution for mixed-stock analysis of 1372 captured fish. The baseline dataset provided sufficient resolution for reliable mixture analysis at regional group level, and also for most of the individual rivers stocks. The majority (76-80%) of the total catch originated from Finnish sea trout populations, 6-9% came from Russian and 12-15% from Estonian populations. Nearly all Finnish trout in the catch were of hatchery origin, while the Russian and Estonian trout were mostly of wild origin. The proportion of fish in the Finnish catches that originated from rivers with natural production was at least one fifth (22%, 19-23%). Two different spotting patterns were observed among the captured trout, with a small and sparsely spotted form being markedly more common among individuals of Russian (28%) and Estonian origin (22%) than among fish assigned to a Finnish origin (0.7%).

Genetic and life-history changes associated with fisheries-induced population collapse.

Over the recent years, growing number of studies suggests that intensive size-selective fishing can cause evolutionary changes in life-history traits in the harvested population, which can have drastic negative effects on populations, ecosystems and fisheries. However, most studies to date have overlooked the potential role of immigration of fish with different phenotypes as an alternative plausible mechanism behind observed phenotypic trends. Here, we investigated the evolutionary consequences of intensive fishing simultaneously at phenotypic and molecular level in Eurasian perch (Perca fluviatilis L.) population in the Baltic Sea over a 24-year period. We detected marked changes in size- and age-distributions and increase in juvenile growth rate. We also observed reduction of age at sexual maturity in males that has frequently been considered to support the hypothesis of fisheries-induced evolution. However, combined individual-based life-history and genetic analyses indicated increased immigration of foreign individuals with different life-history patterns as an alternative mechanism behind the observed phenotypic change. This study demonstrates the value of combining genetic and phenotypic analyses and suggests that replacement or breakdown of locally adapted gene complexes may play important role in impeding the recovery of fish populations.

Discovery and application of insertion-deletion (INDEL) polymorphisms for QTL mapping of early life-history traits in Atlantic salmon.

BACKGROUND: For decades, linkage mapping has been one of the most powerful and widely used approaches for elucidating the genetic architecture of phenotypic traits of medical, agricultural and evolutionary importance. However, successful mapping of Mendelian and quantitative phenotypic traits depends critically on the availability of fast and preferably high-throughput genotyping platforms. Several array-based single nucleotide polymorphism (SNP) genotyping platforms have been developed for genetic model organisms during recent years but most of these methods become prohibitively expensive for screening large numbers of individuals. Therefore, inexpensive, simple and flexible genotyping solutions that enable rapid screening of intermediate numbers of loci (approximately 75-300) in hundreds to thousands of individuals are still needed for QTL mapping applications in a broad range of organisms. RESULTS: Here we describe the discovery of and application of insertion-deletion (INDEL) polymorphisms for cost-efficient medium throughput genotyping that enables analysis of >75 loci in a single automated sequencer electrophoresis column with standard laboratory equipment. Genotyping of INDELs requires low start-up costs, includes few standard sample handling steps and is applicable to a broad range of species for which expressed sequence tag (EST) collections are available. As a proof of principle, we generated a partial INDEL linkage map in Atlantic salmon (Salmo salar) and rapidly identified a number of quantitative trait loci (QTLs) affecting early life-history traits that are expected to have important fitness consequences in the natural environment. CONCLUSIONS: The INDEL genotyping enabled fast coarse-mapping of chromosomal regions containing QTL, thus providing an efficient means for characterization of genetic architecture in multiple crosses and large pedigrees. This enables not only the discovery of larger number of QTLs with relatively smaller phenotypic effect but also provides a cost-effective means for evaluation of the frequency of segregating QTLs in outbred populations which is important for further understanding how genetic variation underlying phenotypic traits is maintained in the wild.

Genetic diversity and differentiation of sea trout (Salmo trutta) populations in Lithuanian rivers assessed by microsatellite DNA variation.

The genetic diversity and differentiation of sea trout were studied in three river basins in Lithuania: Akmena-Dane, Bartuva, and Nemunas. A total of 282 individuals were genotyped at eight microsatellite loci. A similar level of genetic diversity was found in all of the populations studied: mean allelic richness ranged from 3.64 to 5.03, and average expected heterozygosity ranged from 0.588 to 0.721. Significant genetic divergence was observed among the different river basins as well as between populations within the drainages. All pairwise F (ST) values were highly significant, ranging from 0.027 to 0.197. The analysis of molecular variance showed rather weak hierarchical population structuring within the Nemunas basin, which may be explained by extensive gene flow among different river basins or, alternatively, reflect the influence of artificial breeding. Information on genetic diversity and differentiation of the Lithuanian sea trout populations will be useful for future management decisions.